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Strong gravity Lense-Thirring Precession in Kerr and Kerr-Taub-NUT 

spacetimes 

Chandrachur Chakraborty* and Parthasarathi Majumdar^ 
Saha Institute of Nuclear Physics, Kolkata 700064, India 

An exact expression derived in the literature for the rate of dragging of inertial frames (Lense- 
Thirring (LT) precession) in a general stationary spacetime, is reviewed. This formula is valid beyond 
the weak field approximation. When used for the Kerr metric, it leads to the LT precession frequency 
in the strong gravity regime appropriate to compact gravitating objects like rotating neutron stars 
and black holes. Numerical values of the precession rate are computed for a few known cases of 
pulsars and compared to the precession rates in the weaker gravity regimes of the earth and the 
sun. We also derive the exact LT precession rates in Kerr-Taub-NUT and Taub-NUT spacetime 
to show that in the case of spherically symmetric, zero angular momentum NUT spacetimes, the 
frame-dragging effect does not vanish. 
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I. INTRODUCTION 

Stationary spacctimcs with angular momentum (rotation) arc known to exhibit an effect cahed Lense-Thirring 
(LT) precession whereby locally inertial frames are dragged along the rotating spacetime, making any test 
gyroscope in such spacetimes precess with a certain frequency called the LT precession frequency [1]. This 
frequency has been shown to decay as the inverse cube of the distance of the test gyroscope from the source 
[2] for large enough distances where curvature effects are small, and known to be proportional to the angular 
momentum of the source. The largest precession frequencies (fi) are thus expected to be seen very close to the 
source (like the surface of a neutron star or a pulsar, or the horizon of a black hole), as well as for spacetimes 
rotating very fast with large angular momenta. 

Earlier analyses of the LT effect [3] assume large distances (r > M, M is the ADM mass of the rotating 
spacetime due to a compact object like a neutron star) for the test gyroscope. Further, the rotating spacetime 
solution is usually approximated as a Schwarzschild spacetime, and the effect of rotation is confined to a 
perturbative term added on to the Schwarzschild metric. Basically, this is tantamount to reducing the spacetime 
metric to a flat metric (Minkowski spacetime) together with the added perturbation. This leads to the standard 
result for LT precession frequency in the weak field approximation, given by [3, 4] 

Qlt = ^[HJ-f)f - J] (1) 

where, r is the unit vector along r direction. In a recent work reported in ref. [5], an alternative approach based 
on solving the geodesic equations of the test gyroscope numerically, once again within the weak gravitational 
field approximation, is used to compute the frame-dragging effect for galactic-centre black holes. 

In another very recent related work [6] , Hackman and Lammerzahl have given an expression of LT precession 
valid up to first order in the Kerr parameter a for a general axially symmetric Plebanski-Demianski spacetime. 
This is obviously a good approximation for slowly-rotating compact objects. The LT precession rate has also 
been derived [7, 8] through solving the geodesic equations for both Kerr and Kcrr-de-Sitter spacetimes at the 
polar orbit in terms of Appell's first hypergeometric function. These results are not applicable for orbits which 
lie in orbital planes other than the polar plane. We understand that observations of precession [9] due to locally 
inertial frame-dragging have so far focused on spacetimes where the curvatures are small enough; e.g., the LT 
precession in the earth's gravitational field which was probed recently by Gravity Probe B. There has been so 
far no attempt to measure LT precession effects due to frame-dragging in strong gravity regimes [10]. 

While most textbook calculations of the LT precession rate focus usually on the weak- field approximation, an 
exception is the classic work of Misner, Thorne and Wheeler, [11]. Here, the orbital angular velocity for locally 
non-rotating observers in a Kerr- Newman spacetime is given in eqn. (33.24) as an exercise. This formula does 
not appear to be restricted to the weak-field approximation. However, from an astrophysical standpoint, it is 
not clear that the computed angular frequency corresponds to what might be measured as the Lense-Thirring 
precession in a strong gravity situation. Because it has been derived in a locally non-rotating frame which 
the authors amply clarify is not a 'Copernican' frame, (A Copernican frame is a local orthonormal tetrad at 
rest (so moving only in the "time" direction determined by the timclikc Killing vector of the spacetime) and 
"locked" to the spatial part of whatever such coordinate system is chosen, so that it is "at rest" with respect 
to the local inertial frames at infinity. Within this frame, an untorqued gyro in a stationary but not static 
spacetime held fixed by a support force applied to its center of mass precesses. Since the Copernican frame does 
not rotate (by construction) relative to the inertial frames at asymptotic infinity ("fixed stars"), the observed 
precession rate in the Copernican frame also gives the precession rate of the gyro relative to the fixed stars.) 
it is not obvious that this angular frequency corresponds to the LT precession rate of a test gyroscope in a 
stationary (but not static) spacetime. A naive limiting procedure does not appear to reduce this frequency to 
the standard weak-field result (1) in 'Copernican' frames quoted in most other textbooks for the LT precession 
rate in a weak gravitational field. For this reason, a rederivation of the LT precession rate in spacetime regions 
of large curvature, from a 'Copernican frame perspective, remains, in our opinion, an important issue. Two 
motivating factors may be cited in support of such a contention. 

First of all, the near-horizon physics of black holes and that of the outer layers of neutron stars emitting 
X-rays from their accretion discs also might need to be reanalyzed in view of the nontrivial LT precession of test 
geodesies in their vicinity. Most extant treatments of accretion disc physics appear to consider these effects as 
unimportant, even though there is not much in the way of a computation to theoretically justify this viewpoint. 
With upcoming X-ray observatories, as well as multi-wavelength strong gravity space probes currently under 
construction, which envisage to make observations of possible frame-dragging effects in strong gravity situations 
in the near future, the need to go beyond the weak field approximation is paramount. A recent work by Stone 
and Loeb [12] has estimated the effect of weak-field LT precession on accreting matter close to compact accreting 
objects. While there arc claims that what has been estimated in this work pertains more to orbital precession, 
rather than precession of a test gyroscope (which remains the classic example of LT precession), it is obvious 
that in the vicinity of the spacetime near the surface of neutron stars (respectively, the horizons of black 



holes) , the large LT precession of test gyroscopes ought to manifest in changes in the predicted X-ray emission 
behaviour originating from modifications in the behaviour of infalling timelike geodesies of accreting matter 
particles due to LT precession. Thus, there is sufficient theoretical motivation to compute LT precession rates 
in the strong gravity regime, in a bid towards a prediction that future probes of the inertial frame dragging 
effect, in such a regime, may correlate with. 

For rapidly-rotating radio-pulsars, even though the so-called 'lighthouse mechanism' provides a first order 
mechanism of pulsating radio emissions, the theoretical framework employed often tends to ignore most general 
relativistic complications, most importantly that of inertial frame dragging near the surface of these neutron 
stars. Without an understanding of LT precession in strong gravity situations, these effects cannot be accurately 
estimated and blithely ignoring them might turn out to be erroneous. 

In this paper, we derive the LT precession rate in a strong gravity regime in a 'Copernican' frame. Our result 
is applicable to all stationary spacetimes, no matter whether they are axisymmetric or not. This result is also 
applicable for all orbits (not only for polar and equatorial orbits), located at various distances and difi^erent 
angles. In this sense, it is rather more general than the result of [11]. Also, the oft-quoted weak-field result (1) 
(in a 'Copernican' frame) for the LT precession rate is readily obtained from this general result, inserting the 
metric for the desired spacetime. 

The paper is organized as follows : in section II, we discuss about some special features of the frame-dragging 
effect in a stationary spacetime within a 'Copernican' frame. In section III, the exact precession rate of the LT 
precession in Kerr spacetime [13], without invoking either the weak gravity approximation or an approximation 
involving the Kerr parameter. The weak-field approximation reported everywhere in extant literature is then 
shown to emerge from our general formulation, in section IV. We then compute in section V, numerical values 
of the precession frequency in the strong gravity regime and compare these to the standard weaker gravity 
numerical values, for various gravitating bodies, like the earth, the sun and several known pulsars, including 
the millisecond Hulse- Taylor binary pulsar. In section VI, we derive the exact LT precession rates in Kerr- 
Taub-NUT and Taub-NUT spacetimes. We conclude in section VII with a summary and an outlook for future 
work. 



II. EXACT LENSE-THIRRING PRECESSION FREQUENCY 

The exact LT frequency of precession of test gyroscopes in strongly curved stationary spacetimes, analyzed 
within a 'Copernican' frame, is expressed as a co-vector given in terms of the timelike Killing vector fields K 
of the stationary spacetime, as (in the notation of ref. [14]) 

or, 

"m = ^Vr'-K^dpK, , (2) 

where, •q^'^'"' represent the components of the volume-form in spacetime and K Sz Cl denote the one-form dual 
to K k, Q., respectively. Note that Cl vanishes if and only if {K A dK) = 0. This happens only for a static 
spacetime. 

Using the coordinate basis form of K = d^, the co- vector components arc easily seen to be K^ = gfiQ. The 
combination in the right-hand-side of eqn.(2) may now be expressed in terms of the metric components 

^ = ;: — ^^^^[.gotj (ffooffw - gokgoi) - goigkigoo,]]dx'^ (3) 



The vector field corresponding to the LT precession co-vector in (3) can be expressed in coordinate basis as 



n 



(a 90l ^ \ 90i a 

goi,j [Ol oo] goO.jOl 

500 / 900 



(4) 



The remarkable feature of the above equation (4) is that it is applicable to any arbitrary stationary spacetime 
(irrespective of whether it is axisymmetric or not); it gives us the exact rate of LT precession in such a 
spacetime. For instance, a Newman-Unti-Tamburino (NUT) spacetime with vanishing ADM mass is known to 
be spherically symmetric, but still has an angular momentum (dual or 'magnetic' mass [15], [16]); we use eqn.(4) 
to compute the LT precession frequency in this case as well. This result is new and rather general, because, 
there is only one constraint on the spacetime : that it must be stationary, which is the only necessary condition 
for LT precession. The utility of this equation is that; if any metric (17^,^) contains all 10 (4 x 4) elements 
non-vanishing, it can be used to calculate the LT precession in that spacetime. In this case, the precession rate 
depends only on non-zero goiii = 1,2,3) components, not on any other non-zero off-diagonal components of 



the metric. So, this equation also reveals that the LT precession rate is completely determined by the metric 
components gofj,{fi = 0,i), and is quite independent of the other components (in co-ordinate basis). 

A careful inspection of (4) reveals that the LT precession rate given by it vanishes under the following 
conditions : 

• If goi vanishes, because this corresponds to a static spacetime. 

• If the line element of a stationary spacetime is of the form, 

ds^ = .goo ix'" ) dx° dx° + 2,gom [x"' ) dx° dx"" + .gy dx' dx^ . (5) 

Here, to = 1 or 2 or 3 and x™ = x^ or a;^ or x^. We can explore it with an example. Suppose, x™ = r and 
then line element would be 

rfs^ = gQQ{r)dx°dx° + 2gor{r)dx°dr + g.,jdx'dx^ (6) 

where, goo{r) , g^rif) rnean that goo,gor £^re the functions of only r. In such a situation, the formula (eqn.4) 
would be 



n = 

But, we see that goo,gor arc the functions of r only. So 



la 9oi ^ \ gor ^ 

gor.j oi do gm,jOi 

goo J goo 



(7) 



^rjWOrJ — ^rjl TT-j — — L) [S) 

for j ^ r and for j = r, e^r/ vanish due to the antisymmetry properties of the Levi-Civita tensor density. 
The 2nd term ^rjigorgooj = £rrkgorgoo,r = for the same reason mentioned in the above. So, this is a very 
interesting situation in where LT precession rate vanishes, even though go; 7^ 0. In the above given example if 
we replace x"^ = 4> instead of x™ = r, the LT precession vanishes. A well-known example is the Kerr metric 
in Boyer-Lindquist coordinates: g^^j, is not a function of (j). It is the function of (500 = go<t>{T,(^)) only r and 9. 
Above all, we know that Kerr metric is axisymmetric. So, it is impossible to bring the Kerr metric into the 
form of the line clement (5), and the LT precession in such a spacetime cannot vanish. 

The general formula (eqn.4) of LT precession rate shows that any stationary non-static spacetime (axisym- 
metric or not) have a non-zero precession rate. One may ask why the stationary spacetimes show non-zero LT 
precession rate even if it is non-axisymmetric. 

A spacetime is said to be stationary if it possesses a timelike Killing vector field ^"; further, a stationary 
spacetime is said to be static if there exists a spacelike hypersurface S which is orthogonal to the orbits of the 
timelike isometry. By Frobenius's theorem of hypersurface orthogonality, we can write for a static spacetime, 

ila'^bL] - (9) 

If ^° ^ everywhere on E, then in a neighbourhood of S, every point will lie on a unique orbit of ^° which 
passes through S. From the explicit form of a static metric, it can be seen that the diffeomorphism defined by 
t — ^ — i (the map which takes each point on each Sj to the point with the same spatial coordinates on S_t), 
is an isometry. The "time translation" symmetry, i ^- t + constant is possessed by all stationary spacetimes. 
Static spacetimes on the other hand, possess an additional symmetry, "time reflection" symmetry over and 
above the "time translation symmetry" . Physically, the fields which are time translationally invariant can fail 
to be time reflection invariant if any type of "rotational motion" is involved, since the time reflection will change 
the direction of rotation. For example, a rotating fluid ball may have a time-independent matter and velocity 
distribution, but is unable to possess a time reflection symmetry [17]. In the case of stationary spacetimes, 
the failure of the hypersurface orthogonality condition (eqn.9) implies that neighbouring orbits of ^^ "twist" 
around each other. The twisting of the orbits of ^° is the cause of that extra precession in stationary non-static 
spacetimes. 

III. EXACT LENSE-THIRRING PRECESSION RATE IN KERR SPACETIME 

One can now use cqn. (4) to calculate the angular momentum of a test gyroscope in a Kerr spacetime to get 
the LT precession in a strong gravitational held. In Boyer-Lindquist coordinates, the Kerr metric is written as, 

ds^ = -\\ — dt^ d(l)dt + i--dr^ + p^dO^ + [r^ + a^ + sm^ ed4)^ (10) 

\ P^ J P^ A V P^ J 



where, a is Kerr parameter, defined as a = -jj, the angular momentum per unit mass and 

p2 = r^ + a^ cos^ 61, A = r^ - 2Mr + a^ (11) 

For the Kerr spacetime, the only nonvanishing g^i = g^^, i = (j) and j, I ~ r,6; substituting these in eq. (4), the 
precession frequency vector is given by 



^LT^ ^ 



.900 \ a I 9o<p 1 ^ 

.900, r .900, r ] Oe — \ .900, e goo.e Cr 

.900 / V ffoo 



(12) 



2v^ 
where, the various metric components can be read off from eqn. (10). Likewise, 

y/^ = p^sme (13) 

In order to make numerical predictions for the LT precession frequency in a strong gravity domain, we need 
to transform the precession frequency formula from the coordinate basis to the orthonormal 'Copernican' basis: 
first note that 

^LT = ^"de + Ve'dr (14) 

nlT = grri^^f + 9ee{^'? (15) 

Next, in the orthonormal 'Copernican' basis at rest in the rotating spacetime, the tetrad vector eo = u, the 
tangent vector along the integral curve of the timelike Killing vector K . In this basis, with our choice of polar 
coordinates, Q.lt can be written as 

VtLT = Vg^^'r + y/g^n'^e (16) 

where, r is the unit vector along the direction 9. For the Kerr metric, 

Q' = -J sine }f,^^f2 , (17) 

p4(p2 _ 2Mr) ^ ' 

f}'' = 2 J cos6> ,, / (18) 

Substituting the values of 51'' and Vl^ in eqn. (16), we get the following expression of LT precession rate in Kerr 
spacetime 



^^T - 2aAfcosg ^3 ;^ f - aMs.ne-JL^__e (19) 



rVA , , ,r . „ p^ -2r^ 

„,— -r — aM sm ^ , „ 

p3 (p2 _ 2Mr) p3 (p2 _ 2Mr) 

The magnitude of this vector is 

«.T(r, 9) = ^3^^;^^^^^^ [4A.^ cos^ 9 + {p^ - Ir^f sin^ 9\ - (20) 

This is the LT precession rate where no weak gravity approximation has been made. It should therefore be 
applicable to any rotating spacetime like pulsars, rotating black hole etc. 



IV. WEAK FIELD LIMIT 

For large distances, the Kerr metric is approximated as a Minkowski metric with the cross term {g^td(j)dt), 
that is 

,9 ,9 4Af a sin^ 6* , , , ,„, , 



and it has been shown that [3] 



^LT = ^[HJ.f)f - J] (22) 



For large distances and in weak gravitational fields (where we can take r > M and r > a), the 2nd term of 
eqn.(19) reduces to 



, .. J sine {-l + ^cos^9)f 



r3(l + ^ cos2 9)1 (1 + f^ cos2 - 2M) 
J sin 9 ^ 



(23) 



Similarly, from the 1st term of eqn.(19) one obtains the following, 



2Jcos6'.r^/l + §^ - ^^ 



r^(l + 2j cos2 6')5 (1 + fj cos2 9-^) 



r5 

2 J COS 9 



(24) 



It follows that 



^LT{r, 9)^ ^ 2 cos 9f + sin 99 (25) 



where, 9 is the colatitudc. The resemblance of this equation with eq. (1) is unmistakcablc. 

The LT precession for a general stationary metric in the weak field limit may also be derived from (4). In 
this approximation, 



^1 X 9ai ^^ ^ 

5oo — —>^,9ij — Oij, — << i, 
500 



Under these conditions eqn.(4) reduces to. 



^ - ■^<^iji9oi,j9i (26) 



As 6; ~ 9/ , a test gyroscope precesses in such a weak gravitational field with an angular velocity 

n^^Vxg, (27) 

where, g= (5017502,503)- Here, 1,2,3 indicate the space components in that spacetime. 

V. PRELIMINARY COMPARISON WITH OBSERVATIONAL DATA 

Let us take, as an example for the LT precession in strong gravitational field, the case of the Hulse- Taylor 
binary pulsar (also called PSR B1913+16) in the constellation of Aquila, located roughly 21,000 light-years far 
from the earth. Its mass, radius and rotational time period (T) are, respectively, 

M = I.UMq = 2120.95 m, 
Rfn 10 km(assuming), 

T = 59.03 X 10-3 s (28) 

respectively. Here, Mq indicates the mass of the sun. Using these numerical values, we can easily calculate 
the angular velocity (w) of the pulsar 

^ 105.96 rad s"^ 



59.03 X lO--'' 

= 3.53 X 10"^ m-i (29) 

and Kerr parameter (taking moment of inertia / w QAMR^, [18]) 

a « 14.13 m (30) 

All numerical values are in gravitational units with G = c = 1 . 

Now, let us imagine a test gyroscope very close to the surface of the pulsar, in order to that we may capture 
the strong gravitational effects for LT precession. So, we take r -^ R. 



Massive objects 


i^ strong LT{r Sid S ^) 


i^weakLT{ra.d S ) 


Earth 


4.02 X 10-"" 


4.02 X IQ-'* 


Sun 


4.56 X 10^"^ 


4.56 X IQ-"^ 


PSR B1257+12 


239 


178 


PSRB1913+16 


24 


18 


PSR J1614-2230 


712 


463 


PSR J1748-2021B 


274 


121 



TABLE I. Comparison between strong gravity and weak gravity LT precession (all these values are calculated at the 
pole {9 = 0) and near the surface (r — >■ R) of that massive object). 



Substituting the above values in eqn. (20) for 9 = (which indicates that we are measuring the LT precession 
at the pole of the pulsar), we get 



n 



strong LT 



2JRc 



(i?2 + a2)|^^2_2Afi?- 



24 rad s 



(31) 



or, 



n 



strong LT 



0.22 



(32) 



On the other hand, if we use our well-known weak field approximation (cqn. no 25) for LT precession rate 
[19, 20], we get 



2Jc 



or, 



18 rad s"^ 



(33) 



0.17 



(34) 



So, we can easily see that there is really a substantial (^ 30 %) difference between the weak and strong 
gravity LT precession rates in this case. But, if this precession rate for a gyroscope or a satellite due to the 
rotation of the earth or the sun, the numerical results from the weak and strong gravity LT precession rate 
formulas appear to agree to better than 1 %. So, for the earth or the sun, it does not matter whether one takes 
the strong gravity eqn. (19) or the weak gravity LT eqn. (25). 

In Table 1, we offer a comparative numerical study between calculated values flstrongLT and flweakLT for a 
few recently discovered pulsars, with the sun and the earth for comparison. The effect of strong gravity is as 
high as above 100 % in one particular example, and certainly above 30 % in general for all pulsars. 

During the ^strongLT calculations, the backreaction effects of the binary companion of a particular pulsar 
can be usually neglected, because, the distance of this companion is in the order of 1 lt-sec(10^ km) or more 
than that and Vl^t is inversely proportional to p^. The strong gravity LT effect computed which we focus 
on is significant near the surface of pulsars and near black hole event horizons. If the binary system under 
consideration has a non-compact companion of the pulsar which is reasonably distant, backreaction effects of 
this companion on LT precession of test gyroscopes may be safely ignored. 

We should be careful in pointing out, however, that the table of comparisons above is merely to convey the 
importance of strong gravity effects in the rate of dragging of inertial frames in the case of compact gravitating 
objects; the actual numbers may have errors, both theoretical and observational, which have not been fully 
analysed yet. 



VI. LENSE-THIRRING PRECESSION IN KERR-TAUB-NUT SPACETIME 

The Kerr-Taub-NUT spacetime is geometrically a stationary, axisymmetric vacuum solution of Einstein 
equation with Kerr parameter (a) and NUT charge (n). If the NUT charge vanishes, the solution reduces to 
the Kerr geometry. The metric of the Kerr-Taub-NUT spacetime is 



ds' 



4(* 

pi 



Ad^f + i^dr^ 



p^de^ + \ sin^ e{adt - Bdcpf 



(35) 



With 

A = r^ - 2Mr + a^ - n^,p^ = r^ + (71-00086*)^, 

A^asm^e^2ncos9,B ^r'^ + a^ + n^. (36) 

As the spacetimc has an intrinsic angular momentum (due to Kerr parameter a), we can expect a non-zero 
frame-dragging effect. We get from eqn. (12), the LT precession rate in Kerr-Taub-NUT spacetime is 



F^KTN _ v ^ 



a sin 9 



P 



M r 



p2 _ 2Mr — "n? p'^ 



(37) 



p2 — 2Mr — V? p'^ 

In contrast to the Kerr spacetime, where the source of the LT precession is the Kerr parameter a, the Kerr- 
Taub-NUT spacetime has an extra somewhat surprising feature : the LT precession does not vanish even for 
vanishing Kerr parameter a = 0, so long as the NUT charge n 7^ 0. This means that though the orbital angular 
momentum (J) of this spacetime vanishes, the spacetime does indeed exhibit an intrinsic spinlike angular 
momentum (at the classical level itself) which we discuss below in more detail. One can show that inertial 
frames are dragged along this orbitally non-rotating spacetime with the precession rate 



nl^ = "^^°=° f (38) 

where, p^ = r^ + 71^. Notice that the precession rate is independent of 9 and also that it vanishes when the 
NUT charge vanishes, as already alluded to above. In fact, for a = n = 0, the Kerr-Taub-NUT spacetime 
reduces to the static Schwarzschild spacetime which of course does not cause any inertial frame dragging. We 
consider this curious phenomenon in somewhat more detail in the next subsection. 

A. Lense-Thirring precession in NUT spacetime 

The Taub-NUT spacetime is geometrically a stationary, spherically symmetric vacuum solution of Einstein 
equation with NUT charge {n). The Einstein- Hilbert action requires no modification to accommodate this 
NUT charge or "dual mass" which is perhaps an intrinsic feature of general relativity, being a gravitational 
analogue of a magnetic monopole in electrodynamics [21]. 

Consider the line element (of NUT spacetime), which is presented by Newman et. al.[15, 22, 23] 



ds^ = -f{r) 
where. 



dt -|-4nsin^ -c 



+ -r^dr^ + (r^ + n^){d9^ + sin^ 9d(l)'^) (39) 

fir) 



r^ -2Mr-n'^ 

/(O = — ZT-zr. — (40) 



Here, M represents the "gravitoelectric mass" or 'mass' and n represents the "gravitomagnetic mass" or 'dual' 
(or 'magnetic') mass of this spacetime. It is obvious that the spacetime (39) is not invariant under time reversal 
t — > — i; signifying that it must have a sort of 'rotational sense', once again analogous to a magnetic monopole 
in electrodynamics. One is thus led to the conclusion that the source of the nonvanishing LT precession is 
this "rotational sense" arising from a nonvanishing NUT charge. Without the NUT charge, the spacetime is 
clearly hypersurface orthogonal and frame-dragging effects vanish, as already mentioned above. This 'dual' 
mass has been investigated in detail in ref. [24, 25], who also refer to it as an 'angular momentum monopole 
[16] in Taub-NUT spacetime. This implies that the inertial frame dragging seen here in such a spacetime can 
be identified as a gravitom,agnetic effect. 

In the Schwarzschild coordinate system, /(r) = at 



r± 



M ± yjKP + n2 (41) 

r± are similar to horizons in this geometry in the sense that /(r) changes sign from positive to negative across 
the horizon and the radial coordinate r changes from spacelike to timelike. But is r = r_|_ an event horizon in 
the sense of the event horizon of Schwarzschild spacetime ? We shall focus on this issue momentarily. For the 
present, we note that the LT precession rate (which can be easily obtained from eqn. (38) also) is given by ^ 

^MTN ^ n{r--2Mr^u-)h^ ^^^^ 

(r^ + n^) 2 



^ If wc take the Taub-NUT metric{eqn.39) in more general form 

ds^ = -f{r) [dt - 2n.(cos 9 + C)d<l>f + -^dr^ + (r^ -|- n'^)(de'^ + sin^ edcf)'^) 

our final result for LT precession rate(eqn.42) in Taub-NUT spacetime remains same. Because, this final result is independent 
of C which is a constant. In the eqn. 39, we take the form of the metricf23l by Misner (C = —1). 
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FIG. 1. Plot of nf/^ vs r forn = 3 m & M = 1 TO 



It is clear that ftfj^^ = on r = ?'-!-, in contrast to the LT precession frequency in the standard Kerr spacetime 
which is maximum closest to the event horizon ! Further, if we plot the magnitude of the precession rate as a 
function of the radial coordinate for r > r+, as obtained from (42), one obtains the profile like FIG. 1. 

Thus, the precession rate is maximum very close to the 'horizon' r = r-|_, but it sharply drops for ?■ — > r-|_, 
most likely becoming ill-defined on the 'horizon'. 

As the metric (39) blows up at r = r±. we should perhaps try a different co-ordinate system where it is 
smooth on the 'horizon'. Following [26], wherein an analytic extension of the metric (39) has been attempted, 
one obtains the transformed metric 



ds'^ ^ {r^ + n^){d0'^ + 



F^ [du] 



dvl 



- {2n/r±){u±dv± 



v±du±) cos 



where. 



{nlr^f{ul-vl) 



cos 



(43) 



u± 



v± 



F' 



r± 



1/2 



r — r± 



r± 



1/2 



r± J \ r± 
In u, V co-ordinate system r could be redefined as 



r± 



2r± 



exp 



cxp 



r 
2^ 



r 

r± 



cosh I - — I 
\2r±) 



" exp (^) sinh ( J-) 



(44) 
(45) 
(46) 



r± 



r± 



'T 



cxp — 
r± / \r± 



(47) 



Recall that locally every spherically symmetric four dimensional spacetime has the structure I2 ® S^ where 
X2 is a two dimensional Lorentzian spacetime. In this Taub-NUT case, the attempted analytic extension 
discussed immediately above leads to a vanishing of the two dimensional Lorentzian metric on the 'horizon' 
r = r_|_, in contrast to the Schwarzschild metric. This might be taken to imply that perhaps the null surface 
r = r+ is not quite an event horizon; rather it is a null surface where ingoing future-directed null geodesies 
appear to terminate, as already noticed in ref.s [6] and also earlier in [27]. So, physical effects on this null 
hypersurface might not be easy to compute, as a result of which the apparent vanishing of the LT precession 
on this hypersurface is to be taken with a pinch of salt. 

The NUT spacetime, for the mass Af = is also well-defined (see, for example, appendix of [16]). We can 
also write down the precession rate only for massless dual mass (NUT charge n can be regarded as dual mass) 
solutions of NUT spacetime. This turns out to be 



'^LT 



n(r'^ — n^)'- 
(r^ -f- n^) 2 



(48) 



At, the points r = ±n, the LT precession vanishes akin to the previous case, but the same caveats apply here as 
well. One may plot the precession frequency as a function of the radial coordinate as earlier : Here, we observe 

that for n = 3, the LT precession is starting for r > 3 and continued to the infinity. Setting /-^ = 0, we 
get that ^]^x i^ maximum at r = ^/2n. In our figure (FIG. 2.), this value is r = 3-\/2 = 4.24 m. Now, we are 
not interested for r < 3. Our formulas are not comfortable in that regions and r < r± is also not well-defined 
for Taub-NUT spacetimes. From our precession rate formulas (38,42,48) at dual mass spacetimes we can see 
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that the precession rate {^J^t) i^ ^^e same, starting from the polar region to the equatorial plane for a fixed 
distance. ^J^^ depends only on distance (r) of the test gyroscope from the 'dual mass'. 

In summary, we have noted in this subsection several subtleties of computing the LT precession rate on 
and near the 'horizon' of a Taub-NUT spacetime, and our results are consistent with earlier literature where 
geodesic incompleteness on this null hypersurface has been noted. 



VII. SUMMARY AND DISCUSSION 



The analyses presented above has two important features ; (a) the LT precision frequency of a gyroscope in 
a 'Copernican' frame within a Kerr spacetime is computed without any assumption on the angular momentum 
parameter or indeed the curvature of spacetime. The only comparable attempt in the literature is that in ref. 
[11], which however is not the same computation as ours, and the result is not the same in terms of metric 
coefficients, (b) The result derived in eqn. (4) is in fact valid, not just for axisymmetric spacetimes, but also for 
general non-static stationary spacetimes, once again without any assumptions about the curvatures involved. 
This result, we believe is applicable to a very large class of strong gravity systems. 

The substantial difference between the LT precession frequency, arising in strong gravity regime and the 
standard, weak field precession rate for inertial frame dragging ought to provide a strong motivation for their 
measurement in space probes planned for the near future. The fascinating world of gravitational effects asso- 
ciated with strongly gravitating compact objects may provide the best yet dynamical observational signatures 
of general relativity. 

At this point we refer to a recent paper by V. Kagramanova et al. [27] where it is claimed that the Lense- 
Thirring effect in fact vanishes everywhere /every when in the Taub-NUT spacetime. In that paper, timelike 
geodesic equations in this spacetime are investigated. The orbital plane precession frequency (fi^ — VLq) is 
computed, following the earlier work of ref. [28, 29], and a vanishing result ensues. This result is then 
interpreted in ref. [27] as a signature for a null LT precession in the Taub-NUT spacetime. 

We would like to submit that what we have focused on in this paper is quite different from the 'orbital plane 
precession' considered in [27]. Using a 'Copernican' frame, we calculate the precession of a gyroscope which 
is moving in an arbitrary integral curve (not necessarily geodesic). Within this frame, an untorqued gyro in 
a stationary but not static spacetime held fixed by a support force applied to its center of mass, undergoes 
LT precession. Since the Copernican frame does not rotate (by construction) relative to the inertial frames 
at asymptotic infinity ("fixed stars"), the observed precession rate in the Copernican frame also gives the 
precession rate of the gyro relative to the fixed stars. It is thus, more an intrinsic property of the classical spin 
of the spacetime (as an untorqued gyro must necessarily possess), in the sense of a dual mass, rather than an 
orbital plane precession effect for timelike geodesies in a Taub-NUT spacetime. The dual mass is like the Saha 
spin of a magnetic monopole in electrodynamics [21], which may have a vanishing orbital angular momentum, 
but to which a spinning electron must respond in that its wavefunction acquires a geometric phase. 

More specifically, in our case, we consider the gyroscope equation [14] in an arbitrary integral curve 



V„5 =< S.a > u 



(49) 



where, a = V„u is the acceleration, u is the four velocity and S indicates the spacelike classical spin four vector 
5" = (0, S) of the gyroscope. For geodesies a = => V^S" = 0. 

In contrast, Kagramanova et al. [27] consider the behaviour of massive test particles with vanishing spin 
S" = [30], and compute the orbital plane precession rate for such particles, obtaining a vanishing result. We 
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are thus led to conclude that because two different situations are being considered, there is no inconsistency 
between our results and theirs. 

We may reiterate that the preliminary numerical results on LT precession rates for various systems is intended 
to serve as motivation for prospective measurements in strong gravity situations, as also for further theoretical 
work towards understanding the emission mechanism of pulsars and x-ray emission from black holes and neutron 
stars. We expect nontrivial modifications to arise from incorporation of frame-dragging effects in the theoretical 
analyses of these phenomena. We hope to report on this in the near future. 

There arc additional avenues of further work currently being explored : the most general axisymmetric solu- 
tion of Einstein's equation given by the Plebaiiski-Demiahski metric is being investigated for an understanding 
of the LT precession in this case [31]. 
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